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Characterising Body Tissues 



Field of the Invention 

The present invention relates to methods for the characterisation of body tissue. More 
specifically, the invention is concerned with the characterisation of body tissue as normal 
(e.g. healthy) or abnormal (e.g. pathological). The invention has particular, although not 
necessarily exclusive, applicability to the diagnosis and management of cancer, 
including breast cancer. 

Background 

In order to manage suspected or overt breast cancer, tissue is removed from the patient 
in the form of a biopsy specimen and subjected to expert analysis by a histopathologist. 
This information leads to the disease management program for that patient. The 
analysis requires careful preparation of tissue samples that are then analysed by 
microscopy for prognostic parameters such as tumour size, type and grade. An 
important parameter in tissue classification is quantifying the constituent components 
present in the sample. Interpretation of the histology requires expertise that can only be 
learnt over many years based on a qualitative analysis of the tissue sample, which is a 
process prone to intra observer variability. 

Despite the relative value of histopathological analysis, there remains a degree of 
imprecision in predicting tumour behaviour in the individual case. Additional techniques 
have the potential to fine-tune tissue characterisation to a greater degree than that 
currently used and hence will improve the targeted management of patients. 

There remains in particular a need for techniques that can be used to characterise tissue 
to a greater degree in vivo, in order that the need for often painful and distressing 
biopsies can be reduced. 

Summary of the Invention 

In general terms the present invention is based on a recognition that Compton scattering 
densitometry techniques can be used in the analysis of body tissue to very effectively 
discriminate healthy and abnormal or diseased tissue and to discriminate types of 
abnormal tissue. Moreover, Compton scattering has been recognised as having 
potential for application to in vivo tissue characterisation techniques. 



1843335-1 



1 



( 



The invention provides a method for analysing and / or characterising body tissue, the 
method comprising: 

obtaining Compton scatter data measured from a body tissue sample on which a 
penetrating (e.g. X-ray) radiation beam is incident; and 

using the data to provide an analysis and/or characterisation of the tissue 
sample. 

Compton scatter results from an interaction that occurs between a photon and an 
electron. For this interaction the electron is assumed to be unbound and acting as a free 
particle. This assumption can be made if the energy of the incident photon is much 
greater then the binding energy of the atom. Figure 1 illustrates the Compton interaction, 
where E 0 is the energy of the incident photon, E-, is the energy of the scattered photon, 
m 0 c 2 is the rest mass energy of the electron and 9 is the scattering angle of the photon 
and cp is the scattering angle of the electron. T is the kinetic energy imparted to the 
electron. 

The electron taking part in the interaction is assumed to be stationary, i.e. the initial 
energy (E e ) and momentum of the electron equals zero. During the interaction the 
photon imparts some of its energy to the electron. The amount of energy transferred 
determines the angle of the recoil of the electron and the angle of the resultant photon. 

The angle and energy of a Compton scattered particle can be accurately calculated 
using the principle of conservation of energy and momentum. From Figure 1 it can be 
seen that the incident photon has energy E 0 = hv and the scattered photon has energy 

E y = hv 1 . Resolving the energy and momentum into parallel and perpendicular 
components gives the important Compton scatter equation 

F 

1 - f r> \ 



1 + 
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K m 0 c j 
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hence a measure of Compton scatter can be made by detecting the appropriate energy 
photons at a given angle. 

In some instances it may be sufficient for the Compton scatter data to be as simple as a 
count of photons detected at a selected angle/energy in a given time period. In other 
instances, it may be desirable to obtain an absolute measure of electron density (or 
some other derived measurement). Particularly in the latter case, the Compton scatter 
data is preferably corrected for attenuation in the tissue sample. 
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One way to compensate for attenuation effects is to use two radiation sources and two 
detectors. This is an approach commonly used in bone densitometry, but is less 
preferable when examining tissue samples, particularly in vivo, because it results in a 
greater dose of radiation. 

A preferred method to correct for attenuation effects is to obtain data representing a 
measure of the directly transmitted x-ray radiation for each Compton scatter 
measurement. This data can then be used to correct the Compton scatter data for 
attenuation in the tissue sample. 

Especially at low angles (less than 90°), it is also important to be able to distinguish 
Compton scatter measurement from the coherent scatter peak. So, where the 
transmitted radiation is to be used to correct for attenuation it is preferable that the 
energy of the scattered photons detected is as close as possible to that of the 
transmitted radiation. This ensures that the attenuation coefficients are not too different 
for the two measurements. The energy of the incident penetrating radiation beam and 
the angle selected for Compton scatter measurement are chosen such that the Compton 
and coherent scatter peaks can be resolved, whilst minimising the separation (i.e. 
energy) of these peaks. This substantially eliminates self-attenuation effects as it allows 
one to assume that the attenuation coefficients in the sample affecting both peaks are 
substantially the same. 

Preferably the data is used as the input to a predefined calibration model that relates the 
Compton scatter data to one or more tissue characteristics (e.g. normal or abnormal). It 
is particularly preferred that the Compton scatter data is used as an input to a 
multivariate model as described in our co-pending UK patent application GB0328870.1 

Brief Description of the Drawings 

An embodiment of the invention is described below by way of example with reference to 
the accompanying drawings, in which: 

Figure 1 illustrates the cylindrical geometry used as the sample holder for the 
measurement of the electron density; 

Figure 2 illustrates the experimental configuration for the Compton scatter 
measurements; 

Figure 3 shows a scatter spectrum from a malignant breast tissue sample; 
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Figure 4 shows a calibration graph of the calculated linear scatter coefficients against the 
counts measured in the Compton scatter peak for the calibration solutions; 

Figure 5 shows a graph of the differential scatter coefficient from experimental data 
against the calculated electron density for the calibration solutions; 

Figure 6 shows a box plot of the electron density results obtained from the tissue 
samples; and 

Figure 7 shows a graph of the electron density values for each tissue type. 
Description of Embodiments 

The invention is exemplified below with reference to in vitro Compton scatter 
measurements from uniform samples of body tissue. The general technique is, 
however, equally applicable to the analysis of non-uniform tissue samples, including in 
vivo applications. 

The experiment was undertaken in two sections; Compton scatter measurements were 
made on all the samples, followed by transmission measurements. This was done in 
preference to the two measurements being made consecutively for each sample. This 
method was adopted for two reasons; firstly to ensure consistency of set-up between 
samples through minimising the moving of equipment and secondly to save time. 

Material and Methods 



The angle and energy of a Compton scattered particle can be accurately calculated 
using the principle of conservation of energy and momentum. If the incident photon has 
energy E x = h v and the scattered photon has energy E 2 ~ hv\ Resolving the energy 

and momentum into parallel and perpendicular components gives the important 
Compton equation 



Theory 



(1) 
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where /?* 0 c 2 is the rest mass energy of the electron and 8 is the scattered angle. 
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Consider the cylindrical geometry shown in figure 1. A beam of photons is in the 
direction AB with energy E 1 and a detector is placed at a scatter angle 9 to the incident 
beam. The number of scattered photons, S, with energy E 2 reaching the detector is 
given by 



S cc (Fpjexp^- f A ju x {x)dx^Qxp(- £ f* 2 {x)dx 



(2) 



where V is the volume of scattering material, p e is the electron density of the material in 
the scattering volume, (j r is the attenuation coefficient of photons at the incident energy, 
and \j z is the attenuation coefficient of the Compton photons, scattered through angle 9 
with reduced energy E 2 7 If the incident energy Ei and scatter angle 9 are carefully 
chosen it can be assumed that Ei«E 2 and therefore that |J,= p 2 . Using these assumptions 
it follows that 



exp^- (x^x^j exp^- j* ju 2 (x)dx « exp^- {x)dx^ 



0) 



and equation (2) becomes 



(4) 



From the exponential law of attenuation it is found that 



-y- = exp (- J//, (x)dx) 



(5) 



where l 0 is the incident photon intensity and I is the transmitted photon intensity. 
Therefore by obtaining a measure of l 0 , I and S, the electron density can be found from 



(6) 



where S is the scattered count intensity, T=l/I 0 and k is a constant which includes the 
volume term determined using calibration materials with known electron densities. 



Samples 

A sample set of four different tissue types were examined comprising of 5 fibroadenoma 
(benign), 8 invasive ductal carcinomas (malignant), 4 fibrocystic change (non-malignant 
abnormal) and 5 pure adipose (normal) samples. Each sample was examined at two 
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points. The samples were placed into polythene sample vials of 6mm inside diameter 
and 1mm wall thickness. Although the walls of the vial were relatively thick and would 
cause some attenuation of the scattered beam, these containers were chosen because 
they offered a number of important advantages. 

The sides were completely rigid so the samples could be placed into a vial and lightly 
compressed with a stopper without it distorting. This stopper is to remove any air gaps 
and it also minimises tissue movement throughout the experiment. The containers were 
cheap so each sample could have its own holder for the duration of the experiment, 
making it possible to move the sample and reposition it accurately. The samples also 
needed to be symmetrical about a centre of rotation. 

Method 

The Ka characteristic lines produced by a tungsten target x-ray tube were utilized as a 
monoenergetic source to ensure that the Compton scatter peak was detectable. Using 
this method the Compton and coherent scattered peaks from a recorded spectrum can 
be easily resolved and windowed and the bremsstrahlung background subtracted. The 
desired outcome of the experiment was to be able to resolve the Compton and coherent 
scattered peaks, whilst keeping them as close in energy as possible. The detector 
characteristics dictated that the minimum resolvable energy is about 1 keV. 

The experimental set-up is shown in figure 2. The x-ray beam was collimated to 0.5mm 
diameter, both before and after the sample. This was the smallest beam size viable 
whilst maintaining a reasonable flux. The K a line from the tungsten target of the x-ray 
source (E kot2 =57.97 keV) was used. At this energy a scattering angle of 30° gave a peak 
separation of 1 keV between the Compton and coherent scatter peaks. The scattering 
volume comprises of the tissue contained within the intersecting area of the incident and 
scattered beam. For this beam collimation and scattering angle the entire scattering 
volume was contained within the sample, with no air or polythene from the vial included. 
Each sample was measured for a total time of four hours, with the sample being rotated 
throughout the measurement in order to reduce any errors due to the inhomogeneity of 
the tissues. 

Equipment 

The experiments were performed using a Pantak HF160 industrial x-ray tube. An HPGe 
detector was used in order to produce the energy resolution required to resolve the 
Compton and coherent peaks. The energy resolution was measured to be 0.435 keV at 
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59.54 keV (0.73%). The detector was connected via a pre-amp and an amplifier to two 
single channel analysers, one to record the Compton peak and one to record a 
background region. An observed scatter spectrum of a malignant tissue is shown in 
figure 3. 

The two coherent peaks of the K a i and K a2 W lines can be identified and the two smaller 
Compton scatter peaks can be seen. The Ka 2 Compton peak was windowed over an 
area where there was no superposition of the K a2 coherent peak. This windowed area, 
which was used for the scatter measurements, is also shown in figure 3. The 
transmission measurements for each sample were made by placing the detector at zero 
degrees and recording the photon intensity with and without a sample in position in the 
beam. 

System calibration 

As the composition of the tissues being measured is unknown, the electron density 
measurement system needed to be calibrated. This was carried out by measuring 
substances with a known electron density or one that could be calculated. Five 
substances were chosen in order to produce a calibration curve. 

The solutions chosen were water, iso-propanol, and solutions of potassium hydrogen 
phosphate K 2 HP0 4 . Water and propanol were chosen because they are readily 
available, easy to handle and have a known electron density that is close to that of 
biological materials. The concentration of the phosphate solutions could be varied to 
provide solutions with differing electron densities. In order to have values close to that of 
tissue, solutions of 2%, 5% and 10% were used. 

In order to verify the scatter data for the calibration solutions the linear differential 
scattering coefficient can be calculated theoretically as the composition of these 
solutions is known. The linear scattering coefficient is a measure of the probability that a 
photon of incident energy E will be scattered through an angle 0 and is given by 
equation (7): 




S(x) dojx 
M dQ 



(7) 



where — = 

M V 



(8) 
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where M is the molecular mass of the material, p is the mass density; and N A is 
Avogadro's constant. S(x) is the incoherent scattering factor and the differential 
scattering cross section is denoted with KN for the Klein-Nishina cross section. The 
Klein-Nishina differential scattering cross section for the Compton effect is given by 



da 



dCl 



1 



-i3 «- 



I -f-or(l -cos 9) 



1 + cos 9 



1 + 



a 2 (l — cos&Y 
(l + cos 2 9 Jl + a{l - cos e)] 



(9) 



where E is the incident photon energy and 9 is the photon scattering angle, a is the ratio 
of the incident photon energy to the electron rest mass energy given by 

(10) 

and r 0 is the classical electron radius. 

The Klein-Nishina differential scattering cross section is dependent on photon energy 
and the angle of scatter. It was calculated to be 7.177 x 10" 26 cm 2 /electron for 57.97 keV 
photons at a 30° scattering angle. Using this value and tabulated values for S(x) taken 
from Hubbel! et ai (1975) a value for pcompton for each calibration solution was calculated. 
A graph showing the experimental scatter measurements against the scatter coefficient 
values calculated from equation (7) is shown in figure 4. 

The corrected scatter counts are the counts measured in the scatter peak corrected for 
attenuation and are given by 



or = 



E 



m 0 c 



5„ = 



_ [ S mens " B s ] 
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where S CO rr is the counts recorded in the scatter peak corrected for attenuation. S mea s is 
the number of counts in the scatter peak, B s is the background counts in the scatter 
peak, Imeas is the number of counts in the transmitted peak, B T is the number of 



o 



background counts in the transmitted peak, l 0 is the unattenuated count intensity and B 0 
is the background area for these counts. Figure 4 can be used to convert the corrected 
counts measured into differential scatter coefficients for Compton scatter,p s , where 

(12) 

In equation (12) S corr is the corrected scatter counts as described in equation (11), N is 
the systematic experimental error and k is a constant that is found using the calibration 
curve. 

As the composition of the calibration solutions are known the electron densities of these 
solutions can be calculated using the following formula 

where p is the physical density of the material and Z/A is the ratio of atomic number to 
atomic weight for each element of mass fraction U). Z/A values are tabulated and were 
taken from (Attix 1 986). Figure 5 shows the theoretical electron densities calculated from 
equation (13) plotted against the measured scattering coefficients given by equation 
(12). It can be seen that the two quantities correlate well with a gradient equal to the 
Klein-Nishina cross section which is expected, as for high values of x the incoherent 
scattering factors become equal to Z. 

Results 

Figure 6 shows the results of the electron density measurements that were obtained 
from two points on each sample. The median of each tissue type is shown (thick middle 
line). The interquartile range is contained within the box and the whiskers show the total 
range. 

The graph of figure 4 gives a calibration equation to convert the number of counts in the 
corrected scatter peak into the differential linear scatter coefficient p s - The equation 
given by the graph is 
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ju s =1.737xl0- 7 x + 7.919xl0- 3 



(14) 



where x is the corrected counts in the Compton peak. 

These experimental scatter coefficients are then converted into electron densities using 
the Klein-Nishina cross section. This conversion is shown by the trend line in figure 5. 
The average results are shown in figure 7. 

In figure 7 the values of electron density for standard tissue compositions given in ICRU 
report 44 (ICRU 1989) are also displayed. In this report three separate values are given 
for different tissue compositions. The elemental compositions of these six tissues have 
been given in table 1 . It is worth noting that the values quoted in this report are for 
healthy tissues only, as there is no published data for malignant tissue growths. 



Tissue 


|H 

i 

t 

i 
I 


C 


N 


O 


Other 


Adipose #1 


11.2 


51.7 


1.3 


35.5 


0.1 Na, 0.1 S, 0.1 CI 


Adipose #2 


11.4 


59.8 


0.7 


27.8 


0.1 Na, 0.1 S, 0.1 CI 


Adipose #3 


11.6 


68.1 


0.2 


19.8 


0.1 Na, 0.2 S, 0.1 CI 














Glandular 
#1 


10.9 


50.6 


2.3 


35.8 


0.1 Na, 0.1 P, 0.1 S, 0.1 
CI 


Glandular 
#2 


10.6 


33.2 


3 


52.7 


0.1 Na, 0.1 P, 0.2 S, 0.1 
CI 


Glandular 
#3 


10.2 


15.8 


3.7 


69.8 


0.1 Na, 0.1 P, 0.2 S, 0.1 
CI 



Table 1. The elemental compositions (percentage by mass) of adult tissues (ICRU 
Report 44, 1989) 
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It is usually assumed that malignant tissue has approximately the same structure as 
healthy glandular tissue. This is because tumours are usually within fibrous tissue rather 
than growing in purely fatty (adipose) tissue. 

The final results obtained are displayed in table 2. 



Tissue 


Electron density (e/cm 3 ) 


Benign 


(3.330 + 0.140) x 10* 3 


Malignant 

« 


(3.490±0.147)x 10 23 


Adipose 


(3.281 ±0.1 38) x 10 23 


Fibrocystic 




change 


(3.752+0.158) x 10 23 



Table 2. Experimental values obtained for tissue electron densities 



Each individual measurement is subject to statistical variation. The error a is given 

as: a* = J— where x is the mean number of counts if the reading is repeated N times. 

For the scatter readings each measurement was measured for a sufficient time (4 hours) 
to ensure that the error on the counts was sufficiently low (<0.5%). Due to time 
constraints the readings were not repeated. 

The largest error is associated with the subtraction of the background counts. The 
overall error on the background count calculation is 4.2%. This is shown by the error 
bars in figure 4. Other sources of error are the effect of multiple scatter, the error in 
positioning and the error in repositioning the sample for the transmission measurements. 
There is also a widening of the Compton scatter peak caused by the acceptance angle 
of the pre-detector collimator. None of these other errors have been considered as they 
are difficult to quantify and are small compared to the background subtraction error 
outlined above. 
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Discussion and Conclusions 

The results show that there is a detectable difference between the electron density of 
adipose and malignant tissue, to a value of 6.4%. This difference is consistent with the 
values found by using the adipose and glandular tissue values from ICRU report 44. The 
average value for glandular tissue (ICRU44 glandular#2) is 6.2% higher than the 
average adipose value (ICRU44 adipose #2). 

There has been no composition values published for benign (fibroadenoma) or 
fibrocystic tissues. However the measurements made within this study found a difference 
in the electron density of benign and malignant tissues to the value of 5.6% and a 
difference between fibrocystic change and malignant tissue to be 2.3%. 

It is difficult to verify these results using the literature as there is no published data from 
any previous studies using this tissue type. However the high degree of correlation for 
the calibration solutions (figure 4) shows that the system has a reliable accuracy. This 
inspires confidence in the findings that there is a measurable difference between the 
benign and malignant tissues. 

There is a great deal of evidence to suggest that the metabolism and physiology of 
tumour cells differ greatly to that of normal and indeed benign cells. 

Within a benign tumour growth there is often an increase in cell proliferation but the cells 
themselves are relatively normal. However, in a malignant lesion the structure and 
metabolism of the tumour cells and host tissue have a different biochemical structure 
(Gould 1997). This implies that the increase in the electron density of benign tissues 
compare to normal may potentially be due to an increase in cell concentration rather 
than to changes in composition, as seen in malignant tissues. This is consistent with the 
finding that benign tissues ex-vivo have an electron density which is only slightly higher 
than normal tissues and malignant tissues display a much larger difference. 

Dr Otto Warburg first observed in 1930 that cancer cells have a fundamentally different 
energy metabolism than normal cells (Warberg 1930). Since then research has shown 
that tumour cells undergo anaerobic glycolysis, the process where glucose is converted 
to lactic acid through the process of fermentation. This process is extremely inefficient 
compared to normal cell aerobic respiration. 

Glucose consumption rate has been shown to be proportional to histological grade 
(Vaupel et aL 1989) and high grade tumours can absorb about 40 times more glucose in 
order to supply their high energy demands for increased growth. This process is what 
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makes positron emission tomography (PET) imaging so effective at imaging tumours 
using 18F-FDG, an analogue of glucose. It enables PET to distinguish between benign 
and malignant neoplasms with a high degree of accuracy, as benign tissues do not 
exhibit increased glucose consumption (Brock et a/. 1997). 

Anaerobic glycolysis causes a build up of lactic acid to occur within the tissue. The lactic 
acid (CH3-CH(OH)-CO(OH)) which builds up within the tumour has a high electron 
density compared to the host tissue of 8.2 x 10 23 electrons/cm 3 and so could be 
responsible for the increase in electron density that is measured. There is also an 
increase of ketones and glutamine (Vaupel et ai. 1989) which may also increase the 
overall electron density of tumour tissues. Although no direct measurements have been 
made of the composition of benign and malignant tissues, the above suggests that there 
are significant differences in composition. It is difficult to estimate the precise nature of 
the composition changes, given that there are a number of processes occurring in the 
tissue during tumourgenesis. 

The final tissue type that was examined was fibrocystic change. Although this term 
encompasses a range of histological changes, the majority are characterised by tissue 
fibrosis. This is a scarring process whereby the stromal (connective tissue) component of 
the tissue is increased and collagen accumulates. Although increased mature collagen 
may be seen in a few other benign disease processes in the breast, the most 
pronounced increase probably occurs during fibrocystic change. This may account for 
the finding that this tissue classification had a higher electron density than any other type 
of tissue, even malignancy. When examining the tissues exhibiting fibrocystic change it 
is likely that any fluid filled pockets (cysts) will become dispersed during tissue 
preparation leaving only the dense fibrotic tissue under examination. 

Not only is the present invention useful for assessing increased numbers of 
fibroadenomas, invasive ductal carcinomas and FCC tissues, the present invention may 
also be adapted to assess healthy fibrous tissue and further disease processes. 
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